Coherence assisted resonance with sub-lifetime-limited linewidth 
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We demonstrate a novel approach to obtain resonance linewidth below that limited by coherence 
lifetime. Cross correlation between induced intensity modulation of two lasers coupling the target 
resonance exhibits a narrow spectrum. 1/30 of the lifetime-limited width was achieved in a proof-of- 
principle experiment where two ground states are the target resonance levels. Attainable linewidth is 
only limited by laser shot noise in principle. Experimental results agree with an intuitive analytical 
model and numerical calculations qualitatively. This technique can be easily implemented and should 
be applicable to many atomic, molecular and solid state spin systems for spectroscopy, metrology 
and resonance based sensing and imaging. 



Achieving narrow resonance is of long lasting interest 
to both fundamental and applied sciences such as spec- 
troscopy, precision measurement, metrology, and sensing. 
Resonance lines are often broadened, for example, by 
Dopplcr effects, inhomogcneous local fields, power broad- 
ening, etc. Techniques to reduce these effects include 
laser cooling and trapping of atoms, rephasing via pulse 
sequences and Ramsey type pump-probe method pQ. 
Nonetheless, resonance linewidth is ultimately limited by 
the lifetimes of involved atomic states, which determine 
the natural linewidth or more broadly the lifetime- limited 
linewidth. How to go beyond this limit has been a long- 
standing important challenge. Existing efforts can be 
grouped into two categories. In the first category, sub- 
natural optical resonance is enabled by a much narrower 
resonance coupled to the target resonance ■ There, 
the attainable linewidth cannot beat the natural width of 
the narrower resonance. In the second category, a time 
resolved pump-probe method in the same spirit of the 
Ramsey spectroscopy is used to detect the longer-lived 
subset of an atomic ensemble [5H9] . 

We here report a new method which uses continuous- 
wave (CW) lasers to obtain linewidth far below the coher- 
ence lifetime limit. Our technique utilizes coherence be- 
tween the two lower states in a three-level system exhibit- 
ing electromagnetically induced transparency (EIT) [TU] 
or coherent population trapping (CPT) [IT] . The reso- 
nance of interest is between two lower states and the at- 
tainable linewidth is much narrower than the width set by 
the lifetime of the ground state coherence. The narrowing 
factor is only determined by the laser frequency modu- 
lation parameters, and is limited by technical noise and 
eventually laser shot noise. In a proof-of-principle experi- 
ment, we have observed 1/30 of the lifetime-limited width 
by employing frequency modulated CW lasers. Since this 
technique is easy to implement and only requires a Ra- 




FIG. 1: (Color online) Principle of the sub-lifetime-limited 
resonance, shown between |1) and |2) in a three-level A con- 
figuration. S is the averaged one-photon detuning, and A is 
the two-photon detuning, generated by shifting the energy 
levels with a magnetic field via Zeeman shift (as shown) or by 
varying the frequency difference of E\ and E2. (See text for 
details). 



man/CPT process which exists in many systems such as 
atoms, molecules, and solid state spins, we believe it will 
have broad applications. 

The principle of this technique is illustrated in Fig. [T] 
Two metastable ground states |1) and |2) are resonantly 
coupled to an auxiliary fast-decaying level |3) by two op- 
tical fields Ei and Ei under same frequency modulation 
(FM). As the two-photon detuning A is changed from 
zero to nonzero, cross correlation between converted in- 
tensity (amplitude) modulations (AM) [T^] in E\ and 
Ei exhibits a sharp transition from correlation to anti- 
correlation. The physics can be intuitively understood 
as follows. We examine the individual transmissions of 
Ei and E2 as their frequencies are slowly swept together. 
For A = 0, the two transmission spectra overlap and 
their converted AM are correlated. For A ^ 0, the two 



2 



transmission spectra shift to opposite directions due to 
the sharp dispersion of EIT. This effectively induces an 
offset between the two spectra minima much greater than 
A (expected offset if there were no coherence between |1) 
and 1 2)), creates a considerable range around 6 — where 
the two spectra possess opposite slopes, and leads to anti- 
correlations. Therefore, it is precisely the ground-state 
coherence that enables the sharp transition from correla- 
tion to anti-correlation. This picture also suggests that 
the transition linewidth decreases with reduced FM, as 
shall be verified in our theory and experiment presented 
below. We note that such a transition was previously ad- 
dressed in EIT with noisy diode lasers [T3"lfTS] but with- 
out identification of the sub-lifetime-limited width tran- 
sition characteristic. 

The dynamics of the generic three-level A system is 
described by the Hamiltonian: 



i? = n 1 |3>(l|+fi 2 |3)(2|-J|3)<3| + -(|l)<l 



\2}{2\)+H.c, 
(1) 

where fii,2 are Rabi frequencies of E± and E2, and 
H.c. represents the Hermitian conjugate. Here, f2x,2 — 
D, ri 2 e lX sm vt with v and A respectively the modulation 
frequency and modulation depth. The decay rate of |3) 
is r (3> A), and the transverse and longitudinal relax- 
ation rates between ground states are 72 and 71 . The zero 
time lag intensity cross correlation between E\ and E 2 
is defined by «?( 2 )(0) = (Sh(t)SI 2 (t))/y/{(Shf){(SI 2 ) 2 ), 
where 8I12 denote intensity fluctuations of the two trans- 
mitted fields after atoms. The positivity (negativity) of 
g( 2 ) (0) means that two output fields are correlated (anti- 
correlated). For an optically thin medium, </ 2 )(0) can be 
expressed as 5 ( 2 )(0) = (6 p^S p| 2 )) / V \(Sp 31 ) 2 ) ((5 p 32 ) 2 ) 
[13] , where (•) is a time average, and 5p l 31 and 8p l 32 are 
the imaginary parts of fluctuations away from the steady 
states of the atomic density matrix elements. 

To gain physical insights, we pursued analytical solu- 
tions in both adiabatic and nonadiabatic regimes by as- 
suming £l r \, r 2 — Q and equal decay rates from |3) to |1) 
and 1 2). Without these assumptions, our numerical sim- 
ulations also confirmed the subnatural linewidth charac- 
teristic. When the modulation frequency v is much lower 
than the optical pumping rate and the modulation am- 
plitude \v is much smaller than T, atoms adiabatically 
follow the steady-state solution: 



the first (second) two terms in the square brackets cor- 
respond to the correlated (anti-correlated) intensity. A 
static offset between the minima of the transmission spec- 
tra of Ei and E 2 arises from the anti-correlation terms. 
For 5o <C 1, by Taylor expansion and keeping the first 
and second orders of Sq m Eq. (2), we find that the cor- 
relation (anti-correlation) terms only contain even (odd) 
orders of Sq. This is because the two transmission spectra 
have quadratic and parallel FM-AM slopes for A = 0, 
and have linear but opposite slopes for A ^ 0. Due 
to the p\ 2 factor, the anti-correlation terms are zero at 
A = and then grow rapidly with increasing |A|. When 
|A| = r s 2, anti-correlation and correlation terms cancel 
and g( 2 > (0) = 0. Here T g 2 is the half width at half maxi- 
mum (HWHM) of (0) which takes a simple form when 
7i = 72 and 72 < 2r p : 
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with r° = Equation (3) clearly shows that g (2 \0) 
has subnatural linewidth which increases with increasing 
modulation depth and frequency. 

In the nonadiabatic regime, modulation frequency can 
be much larger than the optical pumping rate and even 
r. In this case, ground-state coherence cannot follow the 
laser FM, but the slopes of transmission spectra of E\ 
and E2 around Sq = remain parallel for A = and 
opposite for A^0. We can obtain the first- and second- 
order responses of the atomic coherence using perturba- 
tion theory when the modulation depth A is small: 
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where pj" are the amplitudes of the nth-order responses. 
Similar to the adiabatic case, even (odd)-order terms cor- 
respond to correlation (anti-correlation) . The HWHM of 
p( 2 )(0) can be derived as: 
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where Sq = S/(T/2) = Xvcos vt/(T/2) is the normalized 

o 2 

one-photon detuning, T p = 2 r(i+<5 2 ) * s ^ e recm ced opti- 
cal pumping rate, p\ 2 and p\ 2 are the imaginary and real 
parts of the ground-state coherence P12 = — ,„ r p , . 
The physics becomes now clear. As seen from Eq. (2), 



where C = j^y J j^+^0j^ > an< ^ ^ * s tne Bessel func- 
tions of the first kind. The linewidth is subnatural pro- 
vided that C > 1, which always holds for relatively small 
modulation depth. 

We performed a proof-of-principle experiment in 87 Rb 
vapor under the CPT configuration, with ground states 
being Zeeman sublevels of 5 2 Si/ 2 , F — 2 and excited 
states 5 2 Pi/ 2 , F = 1. The two photon detuning A is pro- 
duced by the magnetic field. The linewidth of the CPT 
resonance is equal to the ground state decoherence rate 
plus power broadening. We used a narrow- linewidth (< I 
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FIG. 2: (Color online) Measured offset between transmission 
minima of E\ and E2 vs. A for two input laser powers. Lines 
are to guide the eye. 
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FIG. 3: (Color online) AC transmitted signals of the two CPT 
fields (blue and red) for (a) A = and (b) A = 15 kHz. (c) 
An example of g (2) (0) spectrum with FWHM 2.4 kHz, about 
1/30 of the transit width 72/^ = 75 kHz. The input laser 
power was 300 /iW. 



MHz) diode laser, and modulated its frequency via ei- 
ther piezoelectric transducer (PZT) or current. Within a 
three-layer magnetic shield, the vapor cell was placed in- 
side a solenoid which provided a homogeneous magnetic 
field. The cell temperature was maintained at 52°C. The 
laser beam was linearly polarized before entering the cell. 
The left and right circularly polarized components played 
the role of E\ and E2 which were separately detected af- 
ter the cell by amplified photodetectors. The AC signals 
were recorded by an oscilloscope and the cross correlation 
of intensity fluctuations was computed offline. 

We first measured the offset between the minima 
of transmission spectra for E\ and E2. For a fixed 
two-photon detuning A, individual transmissions were 
recorded as laser frequency swept slowly. As expected, 
the two spectra overlap at A = and offset to opposite 
directions for nonzero A (see Fig. [IJ. The offset value 
was then retrieved by fitting the spectra with a double 
Gaussian profile to account for the other excited state. 
The offset v.s. two-photon detuning is plotted in Fig. [2] 
The dispersive shape is in agreement with predictions of 
Eq. ([2]). Also, the offset increases with input laser power 
and then saturates. For 2 mW input power, a maximal 
offset of 60 MHz was obtained at A = 0.1 MHz, corre- 
sponding to an effective amplification factor of 600. 

Next, we looked at the g^(0) v.s. A spectrum in the 
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FIG. 4: (Color online) Dependence of g {2) {0) FWHM on (a) 
input laser power, on (b) the power ratio of the two CPT 
beams for laser power of 300 /iW, and on (c) the modulation 
range for two different laser powers. Full modulation range 
2Xv is 5 MHz for (a) and (b). Lines in (c) are linear fits. In (a) 
and (b), an overall scaling factor was applied to the linewidth 
in the fitting to account for broadening from optical depth 
and other noises. The transit width is 75 kHz. 



adiabatic regime. The laser frequency was modulated by 
alternating the PZT voltage at 10.3 kHz, with the FM 
amplitude within tens of MHz. For A = 0, we observed 
correlated intensity fluctuations which only have the sec- 
ond harmonic of the modulation frequency [Fig.[3^a)]. At 
a small nonzero A, the correlation can become completely 
negative and only the first harmonic appears [Fig. |3][b)]. 
The amplitude of the converted AM is small at zero A 
and much larger at nonzero A because the FM-AM slopes 
are quadratic at A = and linear as A / 0. Figure |3jc) 
gives an example of g^ (0) spectrum with a linewidth of 
2.4 kHz. This is about 1/30 of the lifetime-limited width 
of 75 kHz, obtained by extrapolating the zero-power CPT 
width from the measured width v.s. power curve and 
consistent with the transit broadening estimated from 
the 1/e 2 laser beam diameter of 2.2 mm. 

The dependence of the g^ 2 '(0) linewidth upon the in- 
put laser power and modulation depth was studied in the 
adiabatic regime. As shown in Fig. |4]ja), the linewidth 
is larger at lower laser power, consistent with theoreti- 
cal prediction (dashed line). This is because the offset is 
small for low power and complete anti-correlation cannot 
be reached for the chosen range of laser FM, resulting in 
larger linewidth. Such an abnormal trend is in contrast 
with power broadening in traditional spectroscopy. The 
observed width increase at higher laser power is due to in- 
fluences of the residual amplitude modulation (RAM) in 
the laser. Taking into account of this RAM, our simula- 
tion reproduces the experimental observation [solid line 
in Fig. |4^a)]. The signature of the RAM's influence is 
the appearance of the first harmonic in the correlation 
signal at A = 0, as barely seen in Fig. [3]ja). To show 
the general applicability of this technique, we varied the 
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FIG. 5: (Color online) Measured £ (2) (0) FWHM v.s. modu- 
lation frequency in the nonadiabatic regime. Laser power was 
800 jttW, and the transit width is 72/vr = 63 kHz. 



power ratio of E\ and E2 by placing a quarter-wave plate 
before they entered the cell. We found that the g' 2 '(0) 
width has only a small increase as the field intensities 
become more imbalanced, consistent with the numerical 
simulation [see dashed line in Fig. |4|b)]. This weak de- 
pendence on the ratio of two laser intensities allows to 
optimize the ratio against other detrimental effects such 
as light shifts in precision measuremnts. To check the de- 
pendence of the g( 2 \0) linewidth on modulation depth, 
we tuned the PZT modulation voltage and observed the 
linear dependence of the width (Fig. Efc)), which agrees 
with the trend seen from Eq. |3]). The physics is that a 
larger range of frequency variation "sees" more parallel 
slopes and thus reduces anti-correlation. For higher laser 
power, since the offset is larger (see Fig.||, <?( 2 )(0) width 
is more "immune" to modulation depth increase. 

Finally, we measured the g^ 2 \0) width in the nonadi- 
abatic regime by current modulating a distributed feed- 
back (DFB) Laser with about 2 MHz linewdith. Sub- 
lifetime- limited linewidth was also observed (see Fig. [5]) . 
The dependence of the g^ 2 '(0) linewidth on input laser 
power as well as modulation depth (not shown) follows 
the same trend as in the adiabatic case. For a small mod- 
ulation depth A = 0.2 and modulation frequency from 30 
MHz to 180 MHz, the g (2) (0) width has a small variation 
in consistency with Eq. Q. Width increase at higher fre- 
quency is due to a combination of increasing influence of 
RAM (for signal is weaker at higher frequency according 
to Eq. Q) and slightly larger RAM as we found in the 
laser at higher modulation frequency. Calculation taking 
into account these facts and an overall scaling factor on 
the linewidth to include broadening from optical depth 
and other noises fits the data well. Ability to work at 
both low and high FM frequency is desirable for practi- 
cal applications. For example, in CPT clock applications, 
the laser FM frequency can be chosen to be much higher 
than RF modulation, which makes the two modulations 
compatible. 

The smallest attainable g^ (0) width in practice is lim- 
ited by noises. When modulation depth is so small that 
converted AM signal becomes comparable to noises such 
as laser intensity noise, RAM, electronics noise and even- 
tually laser shot noise, the width starts to increase. In 



our transit-broadened system, we have achieved a min- 
imal linewidth of 1.9 kHz with a signal-to-noise ratio 
(SNR) 16 of 62 for an average time of 16 ms; when 
optimized for the ratio of linewidth over SNR, CPT has 
a linewidth of 166 kHz with SNR of 1500. The factor 
of 3.6 improvement in resolving power is mainly due to 
the distinct nature of the g^ 2 \0) observable. As well rec- 
ognized, laser frequency noise induces amplitude noise 
and is a major noise source in precision spectroscopy in- 
cluding CPT-based applications [HJ[l7]. However, the 
g( 2 )(0) measurement is much less sensitive to noise since 
only the noise locked to the laser modulation frequency 
is detected. An enlarged advantage of this method over 
traditional ones on resolving power is expected with im- 
proved FM techniques |18| . In general, narrower reso- 
nance is at the price of a lower SNR, but linewidth con- 
siderably narrower than traditional ones is desirable be- 
cause unknown systematics in resonance lineshapes make 
high ratio line splitting formidable [5]. Therefore, large 
line narrowing factor combined with resolving power en- 
hancement is the merit of this technique. In addition, 
we have numerically examined the possibility to sharpen 
multi-peaks, and find when CPT fails to completely sep- 
arate two closely spaced peaks at its resolution limit, this 
method can clearly resolve them. 

Creating quantum correlations between two bright op- 
tical fields and also among atomic spins using CPT is 
currently of great interest [HI [5D] . Our work paves the 
way for the realization of such proposals, since inevitable 
laser frequency noise induces classical correlations in the 
same way as studied here |21j , and understanding its be- 
havior is necessary for achieving quantum correlations. 
Furthermore, the sharp transition from correlation to an- 
ticorrelation observed here has a close resemblance to the 
behavior of quantum correlations predicted in [55], and 
are both resulted from the coherence induced sharp dis- 
persion. 

In conclusion, we have demonstrated a conceptually 
new way to obtain resonance linewidth far below the 
lifetime limit. Our methodology might also be suitable 
for other level structures provided that an observable is 
identified which is only sensitive to phase information 
near a resonance center. Narrower linewidth combined 
with improved frequency sensitivity is useful for spec- 
troscopy, precision measurements such as frequency stan- 
dards [231 [53] , magnetometry [T7] and sensing based on 
resonance center location. Furthermore, the ability to 
sharpen multi-peaks enables its application in resonance- 
based imaging |25| . where resonance with narrow width 
and high contrast leads to better resolution. These ap- 
plications should be implementable in a wide variety of 
physical systems such as atoms [23] , molecules [26] , quan- 
tum dots [J7] , diamonds [35] , and doped crystals [JHl 130] 
etc., where Raman and CPT process exist. 
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